State permutations from manipulation of near level-crossings. 
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We discuss some systematic methods for implementing state manipulations in systems formally 
similar to chains of a few spins with nearest-neighbor interactions, arranged such that there are 
strong and weak scales of coupling links. States are permuted by means of bias potentials applied 
to a few selected sites. This generic structure is then related to an atoms-in-a-cavity model that has 
been proposed in the literature as a way of achieving a decoherence free subspace. A new method 
using adiabatically varying laser detuning to implement a CNOT gate in this model is proposed. 
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I. INTRODUCTION 

Two central concerns of quantum information theory 
are representation of data through occupancy of quan- 
tum states and the manipulation of these data. Roughly 
speaking, for these purposes, we desire systems that have 
the two properties: 

1) When the system is left alone it should have some 
eigenstates that are quite localized, in the sense that the 
occupancy of such a state identifies a property of a lo- 
calized physical element (e.g., a particular spin being up 
or down, or a state of a single atom being in one energy 
level or the other. ) 

2) There is a simple mechanics for applying a signal 
from the outside that will move the system from one of 
these "localized" states to another. 

In this paper we discuss two models that show some 
interesting potentialities for this program. The first is a 
novel and rather idealized "spin chain" type of model; the 
second is an "atoms-in-a cavity" model, still quite ideal- 
ized, that has been discussed at some length in several 
previous works by other authors. 

All of the models that we shall consider (and some oth- 
ers in the literature) have a unifying feature that can be 
stated as a nearly-general property of a class of Hamilto- 
nians. We define this class as follows: 

a) We picture the states as points distributed in a 
space, and connect each pair of points that correspond to 
a non- vanishing matrix element of H. We consider only 
cases in which each point is connected to a few neighbors. 
A natural extension that will be required for the atoms- 
in-a-cavity discussion is to include the "environment" as 
one of the points in the space, with dissipative links to it 
from some of the other states. 
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b) Now let there be two scales of strength in these 
connections, strong (S) and weak (W). If there are small 
neighborhoods, or blocs, in this space in which all con- 
nections, as well as the diagonal elements, are weak, 
then there will generally be energy eigenfunctions that 
are nearly confined to these blocs. This can be seen by 
treating the weak links to the outside perturbatively after 
diagonalizing the weak bloc in the absence of connections 
to the outside. Since the denominators in this perturba- 
tion are generally of strong order and the numerators of 
weak order, the perturbation will make little difference 
to the state. We can characterize the perturbation of the 
state as being of order W/S. Likewise, we note that al- 
though other weak blocs somewhere else in the space of 
states will indeed have energies of weak scale, we are in 
general protected from a small denominator by the fact 
that we must go through a weak link to get from the 
first bloc into an intervening strong bloc, and then an- 
other weak link to get to the second weak bloc, where 
we would find energies of the weak order characterizing 
the first bloc. Thus the largest term with a dangerous 
denominator would generally be of order W 2 /SW. (The 
dissipative connections to the environment in the atoms- 
in-a-cavity model to be considered later will come to be 
classified as, effectively, a set of additional strong links.) 

The above statement is, to be sure, neither a theorem 
nor an unfamiliar kind of conclusion. In much more so- 
phisticated work on infinite systems much is known about 
the relations between disorder and localization. But we 
will show that the result above, tempered as it is with 
"generally" (meaning "most of the time" ) is a very use- 
ful design criterion for systems that meet the objectives 
enunciated at the beginning of this paper. 

We will begin by considering cases where the strong 
bloc consists of two states connected by a single strong 
link, such that the Hamiltonian is of the form 
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and the eigenvalues of the two energy eigensates are ±S. 
Adding weak scale connections to other states, will make 
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only a small perturbation to these eigenstates. We may 
generalise this to chains of an arbitrary even number of 
states connected by strong links of similar size - all the 
eigenvalues will be of order strong. However, if we have 
an odd number of states, there will always be a zero eigen- 
value, for example, three states connected via a Hamil- 
tonian of the form 

S 12 \ 

Si2 S 23 , (2) 
S 23 J 

will have eigenvalues {0, ±yJS\ 2 + £23}- Due to the pre- 
sense of a zero eigenvalue, adding weak scale links to 
states outside of this bloc will greatly perturb the eigen- 
states. We regain of a set of eigenvalues, all of strong 
order, however, if we introduce either a coupling Si 3 be- 
tween the 1st and the 3rd states, or diagonal elements of 
strong scale. Dissipative connections, for example, would 
enter as imaginary diagonal terms. 

Once we have reason to take a weak bloc as isolated, 
we then wish to address the movements of probability 
through the action of time dependent external interven- 
tions. In the first chain-of-spins model that we shall 
present, with five sites, this will take the form of a bias 
signal to be applied to the middle site, which is changed 
slowly enough for adiabaticity of a kind to be realized. 
The actual state transformations that we discuss can be 
regarded as fairly elaborate "avoided-crossing" phenom- 
ena. In the atoms-in-a-cavity models of Refs. 0], Q, and 
0| that we consider, after noting how the creation therein 
of a decoherence-free subspace can be qualitatively un- 
derstood using the strong-weak picture, we demonstrate 
an adiabatic method of manipulation by tuning and de- 
tuning that is different from the pulsed, and exactly- 
timed, laser techniques that are suggested in these ref- 
erences. 

The remainder of this paper is structured thus: The 
next section analyses spin-chain models with a W and 
S hierarchy of nearest-neighbor couplings and time- 
dependent biases on selected sites. Section HTD then re- 
examines atoms-in-a-cavity models, while Sec. IIVI pro- 
vides further discussion and a conclusion. 



II. SPIN-CHAIN MODELS 

Our first example uses a chain of five spins. We take 
the spins to be numbered consecutively from left to right, 
with pure exchange interactions built of the operators 
between adjacent spins i and j, 

^'l^W^^W.M (3) 
We choose the Hamiltonian 

H = gi h^+\h^+Xh^+g 2 h^ + f{t)[o-& _ 1], 

(4) 



where we have added a single time-dependent term in- 
volving the operator of the middle spin . We refer 
to the function f(t) as the bias on site #3. Since this 
Hamiltonian commutes with the operator J^. ai , we can 
operate within the set of five states with four of the spins 
up and one spin down. The S-W structure comes from 
taking A <C gi,g 2 , and g\ ^ g 2 . When f(t) = 0, this 
choice has the effect of creating eigenstates of H that are 
very nearly the following, 



1>l = 


(Tl + IT)(TTT)/V2, 


1p2 = 


(U - IT)(TTT)/V2, 


1p3 = 


(run), 


Ipi = 


(TTTXU + ID/A 


^5 = 


(TTTXU - it)/V2. 



We understand these combinations qualitatively by not- 
ing that they would be the exact / = eigenstates if the 
small coupling A were set to zero, and that, treating A 
as a perturbation, the link between the pairs of states 
(#1,#2) and (#4,#5) is of order A 2 . 

Now we ask the question: Of the 120 permutations on 
this set of five states, how many can we implement by 
applying a small series of simple pulses in f(t)7 In form- 
ing these pulses we shall tailor f(t) to the permutation 
that is sought, but we shall insist that f(t) begins at and 
ends at the value zero. By a permutation we mean just 
the reshuffling of the states, modulo phase. This demand 
effectively rules out setting relative phases, which would 
require fine-tuning in any system that is to be considered 
over a period of time. The answer to the question is "all 
120". 

The method uses adiabatic avoided-level-crossing dy- 
namics, with an additional feature that can be embodied 
in the above Hamiltonian, namely that the bias on the 
weakly coupled site can be changed suddenly, without 
affecting the state of the system appreciably, at all times 
when the system is far from any (near) level-crossing. As 
an example, we choose the parameters g± = 30, g 2 = 60, 
A = 1, and begin by defining two bias operations fit), 

f a {t) = {t- t ) 0[t - t O }0[t O +T-t}^ U a (t ), 

hit) = -/„(*) -> U b (t ). (6) 

The signal begins at t = to, grows linearly until t = t + r 
when it is switched off abruptly. We take r = 20 in 
the applications that follow. Any bias f a (t) defines a 
transformation, 

U a (to)=T[exp(-i r + dt'H(t'))]. (7) 

Jto 

We will show that the operation U a effects the per- 
mutation (13), while Ub effects (23). The turn-off of the 
bias at t = t + r leaves the state at time t = t a + r 
evolving with the Hamiltonian of Eq. Q with / = 0, 
so that we are prepared to apply another signal to get a 
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FIG. 1: Level diagrams for the three states that are appre- 
ciably affected by the pulse chain of Eq. ©. The numbers 
on the curves indicate possible initial states, enumerated as 
in Eq. 0. In the adiabatic limit, the system, beginning in 
one of these states, follows the continuous path connected 
thereto. For example, if the system is initially prepared in 
state #1, then it first adiabatically transforms to state #3, 
remains there during the sudden bias jump at t — 40, then 
adiabatically transforms to state #2. States #4 and #5 re- 
main virtually unmixed and at their original energies. The 
energy is is units of A and the time in units of A -1 . 
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FIG. 2: The evolution of probabilities for the Hamiltonian of 
Eq. @ with f{t) from Eq. ©, and the initial state #1. The 
light solid curve is the probability of the system in state #1 
(displaced slightly upward for clarity), the heavy solid curve 
is for #3, and the dashed curve is for #2 (displaced slightly 
downward). The time is given in units of A -1 . 



compound of the permutations. Thus we can write, 
Ux^ 2 = [U a (2r)][U b (T)][U a (0)], 



(8) 



where the time arguments emphasize that the transfor- 
mations are to be executed serially over a total time 
of 3t. Using the composition property of permutations 
(13) (23) (13) = (12), we see that just interchanges 

states #1 and #2, leaving #3 alone. The basis for antic- 
ipating the actions of both the individual and compound 
transformations can be shown in an averted level-crossing 
diagram based on Eq. Q with the bias function given by, 



f(t) = f a (t) + f b (t-T) + f a (t-2r). 



(9) 



Figure 1 shows plots of the energy levels for the first three 
states as they evolve under the above successive changes 
of the bias. 

To show that the states really follow these paths, we 
construct the matrix U(t) — T[exp(— i J Q dt' H(t'))] by 
directly solving the Schrodinger equation with the time- 
dependent bias of Eq. ©. In Figs. 2 and 3, we plot 
|(l,2,3|[/(t)|l)| 2 and | (1, 2, 3|E/(i)|2) | 2 against time in 
this solution, giving the expected behavior, |J7(3r)| 2 = 
I U\^2 1 2 with negligible contamination either from non- 
adiabaticity in the region of small level separation, or 
from transitions induced in any of the three sudden 
changes. The plots show clearly the effects of each of 
the constituent transitions in turn, giving rise to the fac- 
torization indicated in Eq. ©. 

To implement the transformations involving states #4 
and #5, we introduce four more bias functions, 

f c (t) = 2(T-t)0(r-t)0(i)-»E/ c -><352), 




FIG. 3: Same as Fig. 2, but with the initial state taken as 
#2. Together, Figs. 2 and 3 show the interchange operation 
effected by the particular pulses of Q. The time is given in 
units of A -1 . 



fd(t) = -fS) - U d -> (341), 
f e (t) = -2(t - t/2)0(r - t)9(t) 

ff(t) = -fe{t)->U f ->(%>), 



(34), 



(10) 



where, for simplicity, we have chosen t$ = Q. By straight- 
forward multiplication of the operations in Eqs. JfjJ and 
CH, we obtain U d U a -^(U), U c U b -^{25), U a U f U a ^(15), 
U b U e U b -^ (24), U e UfU e -^(Ah). These five composite op- 
erations, together with U a U b U a — >(12) and the four single 
pulse operations U a ,U b ,U e ,Uf, give all of the simple in- 
terchanges. All permutations can be built from these 
interchanges, although in most cases it would be more 
efficient to draw on the cycle-of-three permutation op- 



4 



erators U c and U4 as well, or on further primary pulse 
shape variants that directly embody other operators. 

Now we consider longer chains with nearest neighbor 
couplings, 



H 



ftj 

i=j+i 



(11) 



where <ftj are coupling coefficients, and Eq is the en- 
ergy that gives (Ho) = for the case of all spins up. 
We again adopt two scales of coupling strength: S and 
W. As an example, we consider the pattern of couplings 
{51,2, 52,3, 53,4,...} to be {S, W, W, S, W, W, S, W, 
W, . . . }. Under this scheme, the states #3, #6, #9... are 
only weakly coupled to their neighbors. 

As before, we take only states with one spin down and 
the remainder up. With the W coupling constants turned 
off, the eigenvalues of the states in which the down spin 
occupies one of the large blocs, i.e., {(#1,#2), (#4, #5), 
(#7, #8) ...}, come in pairs of ±5^+1. We assume that 
the S couplings are sufficiently irregular so that the dif- 
ferences in energies between nearest strong blocs, e.g., 
.91,2—54,5, are still of the strong (S) scale. When we add 
the W couplings, the eigenvectors will remain almost lo- 
calized, as can be seen from a perturbation expansion. 
In general, as argued previously, small energy differences 
between two non-adjacent blocs are not a concern, since 
as we move across n additional weak links we can toler- 
ate energy denominators that are smaller by a factor of 

(gw/gs) 2n « 1. 

Thus for most cases the eigenvectors can be ar- 
ranged in a list {(£1, f, T, T, ■■■), (vi, T, T, T, ■■•), (T,T,I 

,T,T,...), (T,T,T,6,T,-), (T,T,t, %,?,•••)> (T,T,T,T,T 

, 1,T>— ), ■■•}■ Here, the £'s and ?7's stand for the sym- 
metric and anti-symmetric two-component eigenstates of 
the strong blocs with one spin up and the other down. If 
the system begins with all of its amplitude in one bloc, 
the amplitude will stay almost within that bloc under the 
evolution governed by Hq. 

Now, generalizing the earlier five-site example, we can 
ask whether biases placed at every third site, #3, #6, 
#9, which are weakly coupled both to the right and 
to the left, have the capability of moving information 
around the whole system. We have looked at an eight- 
site example, with independently manipulable biases at 
sites #3 and #6. In complete analogy with the five- 
site case, we find it is easy, for example, to move one of 
the (#1,#2) eigenstates successively to the blocs (#3), 
(#4, #5), (#6), (#7, #8), again by numerically solving 
the Schrodinger equation for the system. 

As a qualitative summary of these outcomes, we restate 
and sharpen some of the remarks in the introduction: 

(1) If we have a chain of blocs of states connected, one 
to another, through an intermediary state that is weakly 
coupled to both blocs, the dynamics within the blocs will 
be nearly self-contained. That is, if the probability is lo- 
calized within a particular bloc at a particular time, in 
the form of any superposition of the eigenstates belong- 
ing to that bloc, then the probability will remain in that 
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FIG. 4: In a notation close to that of Ref. 0, the states are 
labeled as |atom #1, atom #2} ca v where "cav" is the excita- 
tion of the cavity mode. The states denoted "s" and "a" (for 
the atomic part) are |s, a) = (|1, 2) ± |2, l))/\/~2. The connec- 
tions are those variously induced by the three laser couplings 
and the cavity mode, as explained in the text. 



bloc. Likewise, a bloc of several states in which the mu- 
tual couplings and the diagonal elements are weak will, 
in general, not admix appreciably with strongly coupled 
states to its left or right. 

(2) By putting controllable biases on a weak connection 
site to bring the energy of an associated level to (nearly) 
coincide with a level in an adjacent strong bloc, shifts of 
probability from bloc to bloc can be implemented in an 
orderly and complete fashion. 



III. ATOMS-IN-A-CAVITY MODELS 

An exact analogue to the above mechanics for moving 
a system from one state to another can be applied to 
some models of atoms-in-a-cavity type that have been the 
subject of a series of recentpapers. As an example, we 
consider the model of Ref. [3], with two identical atoms 
in a cavity, and a dynamics that is effectively confined 
to three states, which form a decoherence free subspace 
(DFS). 

In this model, the two atoms have three levels, {0, 1,2}, 
and are placed in a cavity tuned exactly to the 1—2 level 
spacing, with the atoms separately addressable by weak 
laser fields that drive Rabi oscillations. The cavity mode, 
if excited, is also allowed to escape through conversion 
to photons at a partially transmitting wall. This is the 
only decohering process. The connection scheme for the 
relevant states is depicted in FigQ] 

The excitation of the cavity mode from the states in the 
chain that are coupled thereto are the S links, denoted by 
fat double arrows. The thin double arrows are W links, 
and the double arrows leading to D represent the cavity 
mode leakage to the outside. All W links come from the 
laser interactions, which are three-fold: on atom #2, an 
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exactly tuned 0^2 signal and an exactly tuned 1 <-> 2 
signal; on atom #1, just an exactly tuned 1 <-» 2 signal, 
180° out of phase with the corresponding signal on the 
other atom. 

The structure shown in Fig. 4 has a general property 
in common with the models discussed above; there is a 
bloc within which all connections are weak, 



1.1) 



mo 



|l,0)o, 



(12) 



characterized in Ref. [2] as a DFS. The existence of this 
DFS is explained qualitively in the diagram of Fig0] In 
this model both the cavity mode coupling (thick arrows 
connecting states) and the cavity mode decay (labelled 
D) are of order "strong" . The cavity mode decay is nec- 
essary for obtaining the DFS, since in its absence the 
strongly coupled block by itself 



|l,l> 2 ^|s>i^|2,2) 



(13) 



has a zero eigenvalue. With the dissipative links, D, in- 
cluded through imaginary diagonal terms in the Hamil- 
tonian matrix, all three eigenvalues are of "strong" level, 
and are thus approximately decoupled from the three 
states {|1, l)n, \a)o, |l,0)o} which form a DFS, as empha- 
sized in Ref.01- 

Ref. 0] proceeds from this observation to the construc- 
tion of a CNOT gate in which the transformation of 
states is effected by means of accurately timed laser 
pulses. We now demonstrate the action of an alterna- 
tive mechanism for the controlled transformation of the 
states, following an exact analogue of the protocol used 
in the first example in this paper. In this mechanism 
the applied laser fields are slowly changed in frequency, 
rather than turned on and off. 

To show this, we first write an effective Hamiltonian 
operating within the bloc of Eq. I|12l) , as given in Ref. , 
but allowing the laser frequencies to be detuned by a 
small amount, 



n r 



A A t\ 



|l,0) (a|o-e jABt |a)o(l,l|o]+h.c. (14) 



Here 17 is a Rabi frequency, and the detunings of the 
«-> 2 and 1 <-> 2 lasers are given by Aa and Ab respec- 
tively. Transition probabilities in the indicated basis will 
be unaffected by the transformation (acting in the atomic 
space only), ^'(t) — exp(—iAt)^!(t), where 



2A = (A A + A B )|1,0)(1,0| + (A fl - 
- (A A +A B )|1,1)<1,1|, 



A j4 )|a)(a| 



(15) 



giving the new effective Hamiltonian, 

H' eS = A+ ^[|1, 0) <o|o - |a) Q (l, 1| ] + h.c. (16) 

This is exactly the first example of the present paper, 
with the five states reduced to three by taking only the 
symmetric states in place of the (#1,#2) and (#4, #5) 



complexes. We may thus interchange the states, ex- 
actly as before, by adiabatically varying the detunings. 
For example, by adiabatically changing the combination 
Aa + Ab, while keeping A^ — Ab constant, we can ef- 
fect the interchange of the states |1, l)o and |1, 0)o, thus 
performing a CNOT operation. 

The full set of equations describing the evolution of all 
the states depicted in Figgis given below, where we have 
used the Hamiltonian from Ref.Q and added detuning as 
in EalPol 



c |10)o 

C\U)i 

C|22) 
C|20) 

C|10>! 

C|H> 2 



in , ^ 

-y (C|10) — c |H)o + c |22)oJ 

1 

-- (A A - A B )c la)o , 

in ( 1 

"Y I C |a)o + ^ C |2,0)o + C |»! 

--(A A + A s )c|io) , 

in i , . 
— 2 C \ a )° + 2^ A+ B' ^ 01 

in 
2 

V2gc| s)o - kc\h) x , 

9 (c|22) ~ V^Clxi). 
y ( C |a)o - V^C| 20 )oJ ~9C\s) 1 : 

in ( 1 r \ 

-y I -/= c |io) - v2c| 22)o J - 0C|i O)l 

5 C |20) - KC |10)i; 

V2gc\ s)l -2 KC | 11>2 . (17) 



-C|io> + V2.9C|ii) i; 



In these equations, g is the parameter setting the strong 
scale depicted by the fat double-arrows in Fig. while 
k sets the (strong) scale of the dissipative links D to the 
environment (decay of cavity modes). 

We have solved this set of equation numerically, to 
demonstrate the interchange of states |10)o «-> |H)o- We 
begin with a steady situation in which detuning param- 
eters Aa and Ab are substantially greater than fi, and 
the initial state is (very stably) either 1 1 , 1) or 1 1, 0)o- 
Then the detuning is manipulated, using simultaneous 
slow changes of both laser frequencies in order to inter- 
change these two states, in a process similar to those de- 
scribed in section [H] The evolution is shown graphically 
in FigO By contrast, Ref.0] uses the perfectly tuned 
case Aa,b = 0, and takes n = until a pulse turn-on 
time. The pulse is then turned off at exactly the time for 
the interchange |l,0)o <-* |1, l)o to have occurred under 
the influence of precession alone. In either method, the 
transformation represents a CNOT gate, the states |0, l)o 
and |0, 0)o being frozen due to the S-W effect. Note that 
our implementation does not require accurate timing of 
applied fields. 
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IV. DISCUSSION AND CONCLUSION 



We have shown that in systems with appropriate ar- 
rangements of strong (S) and weak (W) couplings, vari- 
able potentials applied to a relatively small number of 
sites can efficiently effect state permutations for spin 
chains with pure-exchange coupling. Using the same ap- 
proach, we found a new way to implement a CNOT gate 
in an atoms-in-a-cavity model discussed by previous au- 
thors @. 

A number of models similar to that of Ref. can be 
found in the recent literature 0, & IE 01 El E3 • F° r many 
of these cases, the generic S-W paradigm developed in the 
present paper provides a unified basis for understanding 
how isolated subspaces are generated. 

As our final comment we note that the S-W non- 
dissipative links, and decay-induced dissipation - the two 
cases specifically studied here - are not the only means 
by which subspaces can become mutually isolated. For 
instance, Ref. [8j and many subsequent works have dis- 
cussed cases in which rapid incoherent scattering can 
freeze a system in a single state or subspace of states. 
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FIG. 5: Fhe top panel shows the survival probability of the 
initial state |10)o, the center panel \a)o, and the lower panel 
1 1 1) o - Fhe horizontal axis displays time in units of f2 _1 . We 
have taken the parameters re = g — 2000fi, Aa = (10 — 
0.2t)fi and Aa = (10 - 0.2i)n. Fhe CNOF opertation is 
implemented by adiabatically changing the detunings Aa + 
Ab while keeping A a — Ab constant. 
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